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ABSTRACT 

A detailed imaging analysis of the globular cluster (GC) system of the Sombrero galaxy (NGC 
4594) has been accomplished using a six-image mosaic from the Hubble Space Telescope Advanced 
Camera for Surveys. The quality of the data is such that contamination by foreground stars and 
background galaxies is negligible for all but the faintest 5% of the GC luminosity function (GCLF). 
This enables the study of an effectively pure sample of 659 GCs until ~ 2 mags fainter than the 
turnover magnitude, which occurs at M™^'^ = -7.60 ± 0.06 for an assumed m-M= 29.77. Two GC 
metallicity subpopulations are easily distinguishable, with the metal-poor subpopulation exhibiting a 
smaller intrinsic dispersion in color compared to the metal-rich subpopulation. 

Three new discoveries include: (1) A metal-poor GC color-magnitude trend. (2) Confirmation that 
the metal-rich GCs are ~ 17% smaller than the metal-poor ones for small projected galactocentric 
radii (less than 2 arcmin). However, the median half-light radii of the two subpopulations become 
identical at 3 arcmin from the center. This is most easily explained if the size difference is the result 
of projection effects. (3) The brightest (My < —9.0) members of the GC system show a size-magnitude 
upturn where the average GC size increases with increasing luminosity. Evidence is presented that 
supports an intrinsic origin for this feature rather than a being result from accreted dwarf elliptical 
nuclei. In addition, the metal-rich GCs show a shallower positive size-magnitude trend, similar to 
what is found in previous studies of young star clusters. 

Subject headings: galaxies: spiral - galaxies: star clusters - galaxies: individual (M104, NGC 4594) 



1. INTRODUCTION 
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Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained [from the Data Archive] at the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS 5-26555. These observations are associated with pro- 
gram #9714. 
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Extragalactic globular cluster (GC) systems provide 
an observational probe of galaxy formation and evolu- 
tion. GCs may have formed even before galaxies were 
fully assembled, as well as duri ng subsequent star fo rma- 
tion episodes (see references in iBrodie fc Strader"l2 006). 
Since the Milky Way serves as a standard against which 
other GC systems are compared, it is particularly im- 
portant to study GC systems in other spiral galaxies to 
establish whether the Milky Way GC system is typical. 

The study of spiral GC systems is complicated 
by their generally smaller GC populations and un- 
certainties in the internal extinction conditions. 
Nonetheless, data from the Hubble Space Telescope's 
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(HST) Wide-Field Planetary Camera 2 (WFPC2) 
have helped establish the properties of GC systems 
in several spiral galaxies, including detailed stud- 
ies of individual galaxie s iKis slcr-Pa tig et al. I 119991: 
Larsen. Forbes & BrodigJ 1200 It iBarmbv fc Huchral 
2001: Larsen 20 00) as well a,s lim i ted-scal e sur- 
veys llForbes. Brod ic & La rsen I 120011: iLarsenl 120021: 
Goudfrooii et al. I >2003j iChandar. Whitmore fc Leel 
2004|) . While these efforts illustrated the important role 
of space-based telescopes in extragalactic GC studies, 
large GC samples have been difficult to obtain for 
all but the most densely populated systems due to 
the small field of view of the WFPC2. In contrast, 
the HST Advanced Camera for Surveys (ACS) has 
higher sensitivity and doubles the areal coverage. It 
is thus a powerful tool for global studies of GC sys- 
tems, as demonstrated by a number of ACS surveys of 
such systems in early- type galaxies ijPeng et al.l 120061: 
IStrader et afl feOOG: Harris et al. 2001). 

Figure ^ shows NGC 4594, the "Sombrero galaxy" , 
which is a nearby (9.0 ±0.1 Mpc or m-M= 29.77 ±0.03, 
see W3.2\ for derivation), nearly edge-on (i=84°), spiral 
galaxy (type Sa) with a populous GC system. The Som- 
brero is the most massive member of a small galaxy 
group. While its classification as a spiral galaxy is ques- 
tionable (e.g. Rhode & Zcpf 2004, call it an SO), the 
Sombrero does appear to be an intermediate case be- 
tween elliptical and late-type spiral galaxies, therefore a 
unique system to analyze. 

The GC systems of both spiral and elliptical galax- 
ies generally possess at least two GC subpopulations 
chara cterized as metal-poor (b lue) and metal-rich (red) 
GCs ijBrodie fc Straderl 1200^. Due to th e notorious 
metallicity-age degeneracv l|Worthev I I1994D , ages and 
metallicities cannot be independently constrained by 
broad-band colors alone. Ages derived from Keck spec- 
tra of 14 Sombrero GCs, sugge st that both subpop ula- 
tions are ~ 11 — 12 Gyr old l)Larsen et al. |l 2002tl. as 
is the case in the Milky Way GC system. [Held et al.] 
(2003') used VLT spectra to derive mean metallicities of 
[Fc/H]M„e = -1.7 and [Fe/H]red = -0.7 for the Som- 
brero's two GC subpopulations. A ground-based, wide- 
field imaging study (Rgc < 20') of the Sombrero by 
iMoretti et alJ l|2003D showed that the red GCs are dis- 
tributed spherically and suggested their association with 
the ga laxy's bulge rather than its disk. Rhode & Zcof 
l|200l used a 36' X 36' wide-field, ground-based im- 
age to show that the GC surface density drops to zero 
at a projected radius of ^ 19' (~ 50 kpc). They 
estimated a total GC population of 1900 and calcu- 
lated a specific frequency, Sjy = 2.1 ± 0.3, wher e 
Sn =N x10°-^(^^+^^) (Harris fc van den BeTghlli^Sl . 
iLarsen. Forbes fc Brodie . (2001.1 found that Sombrero's 
red GCs are, on average, ~30% smaller in size than its 
blue GCs. 

2. DATA REDUCTION 

The Hubble Heritage Project sponsored a six-pointing, 
ACS mosaic of the Sombrero Galaxy (PI: Keith Noll, 
HST PID 9714). Each pointing consisted of 4 x 675s 
exposures with the F435W (B) filter, 4 x 500s with the 
F555W (V) and 4 x 350s with the F625W (R). Figure [T] 
illustrates the ~ 600" x 400" region covered by the im- 
age mosaic. Each pointing was analysed separately. Raw 



ACS images were prepared for analysis with the stan- 
dard STScI "On-the-Fly" data reprocessing system. This 
data pipeline made use of the latest calibration files (e.g. 
dark current reference file) to prepare raw data and cor- 
rected the significant geometric distortion in the ACS by 
"drizzling" together multiple exposures. The calibrated, 
drizzled data products were downloaded and prepared 
for object detection by removing diffuse galaxy light via 
the subtraction of a median filtered image. To detect 
GC candidates on the backgroimd-subtracted images, 
the FIND task in the D AOPHOTII ap erture photom- 
etry software application ijStetson Ill992t) was used. An 
intentionally low detection threshold (3(t) was used to 
maximize the chances of finding a faint object. Spuri- 
ous detections were effectively removed by matching up 
detected objects across the three bands. 

Aperture photometry was carried out on the original 
images using the PHOT/DAOPHOTII routine with a 5 
pixel radius and individual sky values calculated from a 
10 pixel wide annulus that began 30 pixels from the ob- 
ject center. A larger aperture size was not used because it 
would have contributed more noise from the sky than sig- 
nal from the object. To compensate for the light found 
just outside the 5-pixel aperture, the median aperture 
difference (5—10 pixels) from candidate GCs was sub- 
tracted from all objects. Additional aperture corrections 
were applied according to the published 10 pixel (0.5") to 
nominal infinite values from the ACS photom etric perfor- 
mance and calibration ( Siri anni et ani200,'iD . Although 
ignored in the present study, a systematic error is ex- 
pected on the 10 pixel to infinity corrections because, 
unlike the synthetic stars used to calculate these correc- 
tions, extended sources will show light beyond the 10 
pixel aperture radius. From a few extremely isolated 
GCs in a B-band image, this effect is estimated to be 
slight (- 0.03 mag). 

Galactic foreground extinction corrections were ap- 
plied using the reddening values, E(i? — 1/), which cor- 
responded to the center of each of the six pointings in 
the DIRBE dust maps of'Schlegel. Finkb einer fc Davis I 
(1998). To convert the E{B-V) values to the corre- 
sponding extinction values ( Am) the reddening co nver- 
sion factors from Table 14 of lSirianni et al. I l)2005t) were 
used assuming a G2 spectral energy distribution. Tabled 
summarizes the corrections applied for each filter. Cor- 
rected photometry was then transformed from the ACS 
system to the Johnson BV and Cousins R systems using 
the coefficients and technique described in lSirianni et al. I 
{2005). 

Globular cluster sizes were c alculated using the 
ISHAPE algorithm (|Larsen 1 119991) . This algorithm es- 
timates the half-light radius (R^;) by convolving a point 
spread function (PSF) with an analytic profile, varying 
the FWHM until the best match to the object profile 
is found. At the distance of the Sombrero, typical GCs 
(Rfii = 2.5pc) will show angular radii of ~ 0.06", which 
is slightly larger than the pixel scale of 0.05" pixel"^. 
After testing three different PSF derivation methods (de- 
scribed in the following paragraph) an empirical PSF was 
adopted for estimating object sizes. The final analysis 
made use of mean sizes derived from the three available 
bands. See ! j2.1.2l for details of combining the size infor- 
mation from the three bands, size selection criteria, and 
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Fig. 1. — Sombrero galaxy image mosaic from the Hubble Space Telescope Advanced Camera for Surveys. The cross and circle symbols 
illustrate the locations of globular clusters belonging to metal-poor and metal-rich subpopulations, respectively. Those objects found within 
the projected dust lanes were not included in the analysis to avoid uncertain extinction co nditions and confusion with HII regions. The 
region covered by a HST WFPC2 pointing used in an earlier Sombrero GC system study bv lLarsen. i^'orbes fc Brodie I <20M is shown. 



the estimated size accuracy. 

Model PSFs for the HST cameras are generated with 
the Tiny Tim program, which also supplies a charge dif- 
fusion kernel (DK) that must be convolved with the PSF 
after it is re-binned for comparison with an object pro- 
file. Since the PSF varies across the two ACS CCDs, a 
21 X 12 grid of Tiny Tim PSFs was generated for each 
CCD. ISHAPE was directed to convolve these PSFs with 
a;.King (1962) model (concentration parameter fixed at 
TiItc = 30) and estimate the FWHM of the object pro- 
files on the F435W images. This method yielded no ob- 
jects with the zero half-light radius that would be ex- 
pected for foreground stars, suggesting the object pro- 
files are broadened compared to the Tiny Tim models, 
likely from the "drizzling" process. 

Two methods were explored to resolve this problem. 
The first used an empirical PSF constructed from bright 
isolated stars on the ACS images. The second method 
involved the creation of a broader DK to effectively re- 
distribute the light away from the center of the Tiny 
Tim PSFs, thereby accounting for the drizzle broaden- 
ing. This broader DK was created by uniformly populat- 
ing a blank pre-drizzle image ( "fit" ACS data suffix) with 
a number of isolated pixels with non-zero values. After 
these images were drizzled, the broadened signal from 
the non-zero pixels was averaged and convolved with the 
original DK, producing a broader DK. Figure [5] com- 
pares the estimated R/jz using the original and modified 
DK methods with the empirical-PSF sizes. Compared 
to the empirical values, sizes derived from the broader 
DK are larger by approximately 0.3 pixels. The original 
method produced half-light radii of objects resembling 




Empirical PSF B — band Sizes{pixels) 



Fig. 2. — Comparison of ISHAPE F435W size measurements us- 
ing three different PSFs. In both panels, the sizes calculated with 
an empirical PSF are on the abscissa. These values are compared 
to the same objects whose sizes were estimated with Tiny Tim 
PSFs convolved with a modified DK and the original Tiny Tim 
PSFs on the top and lower panels, respectively. Lines are 1:1 rela- 
tions. Empirical PSFs were adopted for the size analysis because 
foreground stars show the expected zero sizes. 
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GCs that are larger by ~ 0.5 pixels relative to the em- 
pirical sizes. Both the modified DK and original methods 
produced non-zero sizes for those objects that appear to 
be stellar in the empirical sizes. Therefore the empirical 
PSF-derived sizes were adopted since these are the least 
broadened and stars are readily identifiable. 

2.1. Globular Cluster Selection Criteria and 
Contamination 

Selection criteria based upon magnitude, color, size, 
shape, and location were used to construct a catalogue 
of Sombrero GC candidates. A description of the lim- 
its imposed on the photometric properties of all detected 
objects follows in the next subsection. In i)2.1.2l the 
derivation of sizes is described and how they were em- 
ployed in the selection process. Section 12.1.31 addresses 
contamination in the Sombrero GC dataset. 

2.1.1. Photometry Selection 

Objects within the image region containing the dust 
lanes were excluded from the GC candidate catalogue 
(see Figure Pi to avoid HII regions, young massive clus- 
ters, and uncertain extinction conditions. Photometric 
completeness tests were performed by generating artifi- 
cial objects (using ADDSTAR/DAOPHOTII) modelled 
upon an empirical PSF generated from isolated GCs. 
Colors of these artificial objects reflected the mean ob- 
served GC colors. Five thousand artificial objects with 
the same luminosity were randomly distributed across 
the six pointings for each 0.1 magnitude interval between 
18 and 28 mag. The same data reduction procedure as 
described in §2.0 (except size analysis) was carried out 
on these artificial fields and the resulting detection rates 
are shown in Figure El The Rgc of the artificial objects 
were recorded during the tests to obtain a more suitable 
completeness function when analyzing the GC luminos- 
ity function (GCLF) at different galactocentric distances. 
The 50% completeness level from the combined six point- 
ings is V ~ 26.2. 

A faint magnitude limit of V = 24.3 (at the 95% 
completeness level) was chosen to exclude contamina- 
tion that begins to contribute very significantly at fainter 
mag nitude s, as shown in Figure 21 and discussed in Sec- 
tion 12.1.31 This magnitude limit allows access to GCs 
~ 2 mag below the observed turnover magnitude of V 
= 22.17 ± 0.06, or approximately 95% of the expected 
GC population (see ti3.2() . Since the number of objects 
with magnitudes near this limit is so small, essentially 
no selection bias is introduced by imposing a selection 
criteria on a single band to a sample which spans a ^ 0.5 
range in colors (i.e. the number of candidate GCs in- 
creases by a negligible amount if B-band rather than a 
V-band magnitude limit is used). 

Objects were selected as GCs if their colors fell within 
the ranges 0.9 < B-R < 1.7 and 0.5 < B-V < 1.1. As 
shown in Figure ^ these color limits correspond to the 
observed GC colors, with a small "buffer" to allow for 
some photometric scatter among the faintest GCs. To 
convert broad-band colors to metallicities a straight line 
was fitted to the [Fe/H] and color valu es of Milky Way 
GCs with E(B-V) < 0.3 (|Harris Ill99(t) . To supplement 
the metal-rich end, s pectroscopic m etallicities derived 
from Sombrero GCs l)Larsen et al. 1 12002 ) were paired 
with the ACS photometric information and included in 
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Fig. 3. — Completeness values from artificial object detection 
rates. Square, triangle, and cross symbols represent complete- 
ness values from the F555W (V) band images for the entire, outer 
(1.5' < Rgc < 3'), and inner (Rgc < 1.5') image regions, respec- 
tively. 
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Fig. 4. — Color-magnitude diagram illustrating the color and 
magnitude cuts employed. Candidate GCs in the plot are selected 
on shape and size, as well as a visual inspection. Horizontal error 
bars represent the average color photometric errors and vertical 
error bars indicate average magnitude errors for eight magnitude 
intervals. Dotted vertical lines at B— V = 0.5 and B— V = 1.1 
indicate the color selection limits, while the horizontal line at V = 
24.3 represents the faint magnitude limit. The faint limit excludes 
the increasing amount of contamination which begins to appear at 
fainter magnitudes. 



the linear fits. Figure shows this data, which yielded 
the following color-metallicity transformations: [Fe/H] 
« -4.90 + 3.06(B-R), [Fe/H] « -4.45 + 4.39(B-V) and 
[Fe/H] « -5.64 + 9.43(V-R). The GC color selection 
limits approximately correspond to the following metal- 
licity ranges: -2.1 < [Fe/H](B_fl) < -fO.3 and -2.3 < 
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Fig. 5. — Diagram illustrating the B— R color to metallicity trans- 
formation. The transformation is represented by the solid line and 
was derived from a linear fit to Milky Way (M W) GCs with G alac- 
tic reddening values of E(B— V) < 0.3 from the HarrisJ 119961) cat- 
alogue (circles) and the met allicity values derived from Sombrero 
GC spectroscopy (triangles; 'Larsen et al. 2001). For comparison, 
the color- metaUicity transformations from Barmbv ct al. (2002, 
dashed line) and the HBruzual &: Chariot . ^2003, dotted line) 13 
Gyr models are presented. The former is consistent with the MW 
GCs, but deviates from the Sombrero data points in the metal-rich 
GC region, while the later is offset from the metal-poor GCs. 



[Fe/H](B_y) < +0.4. 

Detections were only adopted as GCs if their 
FIND/DAOPHOTII roundness values were between 
—0.45 to 0.45 and their sharpness values were between 
0.55 and 0.85. Roundness and sharpness thresholds were 
based upon the inspection of all objects that resembled 
GC in color and size. Finally, GC candidates were in- 
spected for any obvious contaminants (e.g., HII regions, 
disky galaxies) and for duplicate objects lying on the six- 
image mosaic overlap regions. 

2.1.2. Size Estimates and Selection 

As described in ^ object sizes were determined using 
an empirical PSF and the profile-fitting routine ISHAPE. 
To reduce noise from the fitting process, three size values 
were determined for each object by running ISHAPE on 
each of the B, V and R band images. The final size 
was taken to be the mean of the three size values, which 
first underwent an iterative sigma-clipping algorithm to 
remove outliers and systematics as is described in the 
following paragraph. 

The V band sizes were arbitrarily chosen to be the 
reference system. Using the objects which passed the 
aforementioned selection criteria, systematic differences 
between the bands were estimated by fitting a line to the 
V versus B and V versus R size values. These fits were 
used to remove the systematics by forcing the B and R 
sizes to show a 1:1 trend with the V sizes. Objects falling 
more than Scr from the two 1:1 lines were discarded from 
further analysis. A line was fit to the modified dataset to 




V-band Sizes(pixels) 

Fig. 6. — Comparison of sizes from B-band (top) and R-band 
(bottom) images compared with V-band sizes. Both have under- 
gone an iterative sigma-clipping procedure to remove systematics 
and outliers, as is described in the text. Triangle symbols represent 
the objects that fell outside of the 3cr regions indicated by dotted 
lines. Solid line is a 1:1 line for comparison. 



gain a more accurate understanding of the systematics. 
After applying a correction to the B and R sizes with 
the improved systematic estimate, another 3cr clip was 
applied to each of the final 1:1 lines. For the remaining 
objects the mean of the resulting size values was adopted 
as the final size. 

Figure ini compares the size values of the three bands, 
after the above procedure was applied. The majority of 
the objects that were removed in the sigma clipping al- 
gorithm tend to be fainter objects, as shown in Figure 
[7| A significant difference between the observed object 
sizes in at least one of the bands likely indicates the ob- 
ject is not a typical old stellar system, an idea that is 
supported by the fact that a significant portion of these 
objects deviate from typical GC colors (Figure Q). 

Least-squares RMS values of the final size fits are 0.095 
and 0.082 pixels for the V versus B and V versus R 
sizes, respectively. Thus the sizes are accurate at approx- 
imately the 0.1 pixels (or 0.2 pc) level. To remove stellar 
contaminants from the dataset, a lower-limit of 0.3 pixels 
was enforced to account for any sizes that deviated by 3a 
from a size of zero, as shown in Figure |S1 Objects with 
ISHAPE ellipticities (e = 1 — b/a) greater than 0.5 were 
also discarded as likely background galaxies. 

The completeness tests presented at the beginning of 
this section did not account for the size selections, there- 
fore these completeness functions are not perfect repre- 
sentations of the complete Sombrero GC selection pro- 
cess. The true completeness functions likely differ by a 
small amount, since the objects removed in the sigma 
clipping algorithm tend to be unlike GCs (Figure O and 
the number of stars predicted and found in the mosaic 
make up only a small fraction of the data sample (see 
the following subsection). 
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Fig. 7. — Color-magnitude diagram with the crosses representing 
the final GC dataset and the triangles the objects removed by the 
size sigma-clipping described in the text. The objects tend to be 
faint and show colors unlike typical GCS. 



Fig. 9. — Color-magnitude diagram with crosses representing 
the final Sombrero GCs. The triangles are the objects found on a 
control field. 
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Fig. 8. — Object size versus B— R illustrating the adopted size 
and color selection limits. Data are selected on sharpness, round- 
ness, luminosity, and a visual inspection. Those objects below the 
horizontal dotted line at K^i = 0.3 pixels are likely foreground 
stars. Dotted vertical lines represent B— R color limits of 0.9 and 
1.7. 



2.1.3. Contamination Estimates 

An unknown quantity of contamination from back- 
ground galaxies and foreground stars is present in 
this dataset. While the resolving capabilities of the 
ACS instrument make identifying unwanted object more 
straight-forward, it is nevertheless important to estimate 



how much contamination is missed. Reducing an iden- 
tical set of ACS images taken just beyond the "edge" 
of the Sombrero GC system is the ideal method for es- 
timating the contribution from contamination. Lacking 
such a control field, B and V images (2400 seconds in 
B and 1120s in V) from the Hubble D eep Field North 
of the GOODS survey (HST ID 9583; iGiavalisco 6^111 
l2?ffll were utilized. A single region of the survey was 
reduced in exactly the same manner as the Sombrero 
mosaic and sizes were determined on the B image bands. 
The images were examined and found to have few Som- 
brero GC-like objects when the selection criteria were 
applied, as shown in Figure El Even when one increases 
the observed objects found in the control field by a factor 
of six to account for the six pointings in the ACS mosaic, 
the number of background objects classified as Sombrero 
GCs would be only 24 objects, or ~ 4% of the final Som- 
brero GC sample. This investigation also supports the 
faint Sombrero magnitude limit, which was motivated by 
an obvious increase in the background contamination at 
V ^ 24.3 (see Figure corresponding to the objects 
detected in the control field. 

Since this control field was originally selected to avoid 
stars, it is necessary to estima te the contamination frorn . 
foreground stars. From the iBahcall fc Soneiral l|1980ri 
Galaxy model, approximately 40 stars with the same ap- 
parent magnitude and B— V colors as the Sombrero GCs 
are estimated to be within the ACS mosaic region. After 
applying all photometric selection criteria to the objects 
detected in the Sombrero field, approximately 41 objects 
are classified as foreground stars because they fall below 
the lower size limit of Rhi = 0.3 pixels (i.e. there are 
41 objects within the color limits and below the verti- 
cal line in Figure IHJ. This number is consistent with the 
predicted number of stars, thus the Sombrero GC sam- 
ple likely contains a negligible amount of contamination 
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from foreground stars. 

3. RESULTS AND DISCUSSION 
3.1. Color-Magnitude Distribution 

Figure El shows the B— R, B— V and V— R color- 
magnitude diagrams of the 659 objects that met the GC 
selection criteria. To parametrize the GC color distri- 
bution a Bayesian univariate norma l mixture analytical 
method (^Richards on fc Green Il997j) was used via a soft- 
ware routine named "Nmix"^. Peak colors of B— R — 
1.15 and 1.44 ([Fe/H] = -1.38 and -0.49) were found 
with 53% of the GCs assigned to the blue subpopulation. 
As shown in Table |21 the GC subpopulation peak metal- 
licities and blue proportions are similar to the same val- 
ues f rom the B— V and V— R color distributions. A KMM 
test (jAshman et al.llT994() performed on the B— R colors 
supports a bimodal distribut ion over a unim o dal on e at 
the 99.9% confidence level. iRhode fc Zerfl l|2nf)l de- 
tected a Sombrero bimodal metallicity distribution with 
59% blue GCs and subpopulation peak B— R values con- 
sistent with these results. Since their ground-based study 
extended to the "edge" of the GC system at Rgc — 19' 
(while the ACS mosaic has Rgc ^ 6'), and because blue 
clusters tend to dominate the halo, their larger percent- 
age of blue GCs is not unexpected. In the subsequent 
analysis, blue GCs are taken to be those with B— R col- 
ors < 1.30, and red GCs are those with B-R > 1.30. The 
only exception to this rule are the two brightest GCs with 
B— R > 1.30. It is assumed these are in fact "metal-poor 
GCs" that happen to lie at the extreme of a metallicity- 
mass trend (presented in Section [3. 3() and therefore show 
redder colors compared to typical metal-poor GCs. 

Usi ng a conversion from B — R to V— I from Milky Way 
GCs (jForbes fc Fortd 120011) . the blue and red subpop- 
ulation peaks approxi mately correspond to V— I = 0.93 
and 1.13, respectively. iLarsen. Forbes fc Brodie I l|200lD 
found on V and I-band WFPC2 images, color peaks 
of V-I = 0.96 and 1.21. When compared to the GC 
peak V— I co lor versus galaxy luminosity relations of 
iStrader et aP {2004), the ACS Sombrero blue GCs faU di- 
rectly onto the relation, while the red GC peak is slightly 
bluer than expected. The ACS red peak is likely different 
from the past studies because the red part of the color 
transformation is poorly constrained due to the small 
numbers of red Milky Way GCs. 

To quantify the intrinsic dispersions of the subpopula- 
tions, the sample was limited to V < 22.0 and Nmix was 
used to estimate the Gaussian dispersions: c^"/^ = 0.08 
and crf^;^^ = 0.10. These la color ranges, when con- 
verted to metallicity, are consistent with the dispersions 
derived when the color values were first converted to 
metallicity: 4^^^"^ = 0.23 and al^J^"^ = 0.28. Sup- 
porting this analysis, the same procedure applied to the 
B— V and V— R color distributions yielded similar metal- 
licity dispersions as shown in Table [3 

Other studies have similarly exploited the low photo- 
metric uncertainty of the ACS dat a and estimated the 
intrinsic subpopulation dispers i ons llStrader et al. II2OO6I: 
iHarris et all 120061: iPenEf et alJ \2n()K As found here, 
each of these studies detected a larger color spread in 
the metal-rich subpopulation compared to the metal- 



|http: //www. stats .br is .ac .uk/$\sim$peter/Nmix/ 



poor GCs. However, lacking a standardized color-to- 
metallicity transformation needed for a physical interpre- 
tation, it is difficult to make conclusive statements about 
the intrinsic dispersions of the GC subpop ulations. For 
instance, the present analysis agrees with IHarris et alJ 
([2006^) in that a smaller metallicity dispersion is in- 
ferred in the metal-poor subpopulatio n, but a non-linear 
color-to- metallicity transformation in lPeng et al.l l)2006j) 
yielded a larger metallicity dispersion amon g the metal- 
poor GCs compared to the metal-rich GCs. iPeng et alJ 
(2006) briefly explore the physical implications of the 
subpopulation metallicity dispersions and similarly ex- 
press their reservations about making too many physi- 
cal interpretations until a definitive color-to-metallicity 
transformation is established. 

3.2. Luminosity Function 

Figure ^2 shows the Sombrero V-band GCLF over- 
laid with a Student-t5 (^5) dist ribution fit fr om the 
maximum- likelihood algorithm of lSeckerl (|1992() . This 
algorithm accounts for completeness and photometric er- 
rors. The long exposure time and the resolving power of 
the ACS camera enabled contaminants to be identified 
and eliminated with unprecedented accuracy, resulting 
in an essentially pure sample of Sombrero GCs. By inte- 
grating the best-fitting Sombrero GCLF ^5 distribution, 
it is estimated that only 5% of the GCs are missed from 
the magnitude limit of V < 24.3. Such quality and the 
large sample allows the faint end of the GCLF to be di- 
rectly examined, showing for the first time that the entire 
Sombrero GCLF follows a t5 distribution. 

The location where the GCLF distribution peaks, the 
turnover magnitude (TOM), is a useful quantity to com- 
pare the GCLFs of different galaxies. The apparent 
y-band TOM from a t^ fit, m^^*^ = 22.06 ± 0.06, is 
the equivalent to a Gaussian TOM estimate: 
22.03 ± 0.06. A ^5 distribution is adopted because its 
shape appears to better represent the wings of the GCLF 
(confirmed with a KS test). Comparing the t^ fit to the 
GCLF histogram in Figure ITTl it is apparent that the fit 
may not accurately represent the TOM. Indeed, the peak 
of a smoothed density kernel estimate shown in the same 
Figure suggests a slightly fainter TOM might be more 
appropriate. The peak in density estimate corresponds 
to an apparent magnitude of m™*^ = 22.17, which is 
taken to be the Sombrero GCLF TOM. 

TOM values derived from the two other available bands 
and color conversions from the Sombrero mean GC colors 
are consistent with the TOM values found in the V band. 
The ACS t5 TOM is consistent with a previous WFPC2 
Sombrero turnover estimate of m™*^ — 22.10 ± 0.15 
(.Larsen. Fo rbes & Brodie 2(J^) , but the ^5 dispersion is 
broader than the WFPC2 estimate: af^'^^ = 1.08 ±0.06, 
a^FPC2 ^ 0.81 ± 0.10. The difference between the 
past result is likely a result of small number statis- 
tics and/or an increased amount of contamination at 
fainter magnitudes in the WFPC2 dataset. This is 
demonstrated by a "jump" in the estimated WFPC2 
GCLF dispersion to a^^^^^ = 1.02 ± 0.12 when the 
magnitude limit is increased b y half of a magnitude 
(|Larsen. Forbes fc Brodi e 2001). 

Figure 1121 shows a normalized form of the Sombrero 
GCLF in units of solar luminosity. The TOM corre- 
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Fig. 10. — Color-magnitude diagrams showing the globular clusters of the Sombrero galaxy. Color /metallicity peaks, dispersion values 
and subpopulation proportions for each distributions are presented in Table 1^ Error bars are the same as in Figure |2 As discussed in 
i|3.3l a color-magnitude trend is observed in the blue subpopulation where brighter GCs show, on average, redder colors. This is illustrated 
with solid lines from least-squares fits to peak colors of only the four brightest magnitude bins, which are represented by squares. 



sponds to the bend in the distribution. A power-law 
fit of the form dN{Lv)/dLv = (iLy yields an expo- 
nent of a = —2.2 ± 0.1 for the points brighter than the 
bend. This exponent is consistent with observations of 
you ng star clust ers luminosity functions in spiral galax- 
ies (jLarsenll2002jl . Power-law fits to the Sombrero metal- 
poor and metal-rich subpopulations yields exponents of 
a = —2.1 ±0.1 and a = —1.9 ±0.1, respectively. Fainter 
than the bend, power-law fits to the subpopulations and 
the entire sample are consistent with a flat distribution. 

Observational evidence suggests the GCLF of ellipti- 
cals galaxies and Local Group spi rals (M31 and the Milky 
Way) share a universal TOM ijHarris ]|2002 IRichtlerl 
|2003|) . A composite GCLF from Virgo Cluster dwarf 
ellipticals (dEs) ap pears to have a simila r TOM to mas- 
sive Virgo galaxies ijStrader et al. ll2?)0fil) . Table |21 shows 
distance estimates measured to the Sombrero galaxy 
from a variety of sources. If the midpoint between 
the surface brightness fluctuation (SBF) and the plane- 
tary nebulae luminosity function distance is adopted (i.e. 
m-M= 29.77 ± 0.03), then the ACS absolute TOM is 
y[TOM ^ _7 go ± 0.06. 

After applying an correction of —0.16 mag to account 



for a re-calibration of the iTonrv et all l|2001f) SBF dis- 
tances (as discussed in Jensen et al. 2003) the iRichtler I 
((2003|l "universal" GCLF TOM value is M™^^ = -7.35, 
with a dispersion of 0.25 mag. This is almost three tenths 
of a magnitude fainter than the Sombrero TOM. A slight 
discrepancy such as this may originate in GCLF TOM es- 
timation uncertainties like contamination, completeness 
limits, photomet ric calibration, and varying GCLF fit- 
ting techniques ijRichtler 112003(1 . The distance moduli 
are uncertain as well, therefore the Sombrero TOM is 
likely consistent with the universal value. If instead the 
Sombrero TOM is assumed to correspond exactly to the 
universal TOM, a distance modulus of m— M = 29.52 
(8.0 Mpc) is derived. 

GCLFs of the subpopulations are shown in Figure El 
In contrast to what has been found obse rvationally in 
several galaxies fe.g. iLarsen et al.l I2001D and de mon- 
strated theoretically I* Ashman. Conti fc Z eDf'1995'). the 
blue TOM in the Sombrero is not brighter than the 
red TOM: Ay™5,,, = -0.09 ± 0.12 mag (m™ = 

22.11 ± 0.08, (Juue = 1-08 ± 0.08, m^^/^ = 22.20 ± 0.09 
and ared = 1-16 ± 0.08; TOMs are from the peak of 
a density estimate, while the TOM error and disper- 
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Fig. 11. — Observed V-band GC luminosity function, with 
the best-fitting Student-i5 distribution over-plotted as a thick 
curve (turnover magnitude: m^''*^ = 22.06 it 0.06; dispersion: 
(Tt5 = 1.08 ± 0.05). The thin curve represents the smoothed den- 
sity estimate of the distribution, the peak of which is taken to be 



the Sombrero TOM: m™*^ 



22.17 or M'^OM 



-7.60. From 



artificial field tests it is estimated the 95% completeness level is V 
~ 24.3, thus the ACS mosaic enables the detection of nearly all 
GCs brighter than the faint magnitude limit of V = 24.3. 
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Fig. 12. — GC luminosity function in u nits of solar lu minos- 
ity. A Student-t5 distribution fit from the ISecker I 1199211 algo- 
rithm is represented as a solid line. Dashed lines are power-law fits 
{dN(Lv)/d,Lv = f^Ly) *° regions of the distribution, with the 
corresponding power-law exponents indicated next to each of the 
lines. 



Fig. 13. — GC subp opulation luminosity functions. Plots are the 
same as in Figure lTTl with motal-poor (m^J*^ = 22.11±0.08) and 

metal-rich (m^^'^*^ = 22.20 it 0.09) subpopulations in the bottom 
and top panels, respectively. 



sion values are from the Seeker algorithm). The pre- 
dicted difference for metallicity peaks of [Fe/H]= —1.4 
and [Fe/H]= —0.4, and assuming the two subpopu- 
lations have the same old ages, GC mass functions, 
and stellar initial mass function is AyJ^^^^^^^ = —0.30 
l|Ashman. Conti fc ZeSI 1199,^ . The TOM differences 
from other bands are also marginally smaller than ex- 
pected: AB™*:^^^^ = -0.27 ± 0.12 versus the predicted 
^B^i^e-red = -0.52 and = 0.04 ± 0.12 versus 



the predicted ^Ru^e- 



blue— red 



-0.25. 



The difference between the observed and predicted val- 
ues may indicate that the simplifying assumptions made 
about the GCLFs in [Ashman. Conti & Zepf (199^ do 
not accurately represent the Sombrero subpopulations . 
For example, according to the lBruzual fc Chariot 1 1)20031) 
SSP models if the metal-rich GCs are younger than the 
metal-poor GCs by 1 Gyr then from age effects alone a 
metal-rich GC is expected to be brighter by 0.04 V-band 
magnitudes compared to a metal-poor GC of the same 
mass. 

Within the mosaic region of complete radial coverage 
a TOM difference between the inner (0' < Rgc < 1-5') 
and outer (1.5' < R^c < 3.0') G CLFs is not s tatisti- 
cally significant according to the ISecker I l)1992t) algo- 
rithm: AV^^^it = 0.11 ± 0.13. When the subpop- 
ulations are each divided into the same radial bins, a 
statistically significant TOM difference is not observed 

-0.09 ±0.25 

0.23 ± 0.19 for the blue and red 
GCs, respectively. For each of these radial comparisons, 
the corresponding differences in GCLF dispersion val- 
ues are found to be insignificant. A null result for such 
analysis suggests that various radi al-dependant dynam- 
ical evolution mechanisms (e.g. iFall fc Zhang I l2001t 



among either subpopulation: AV^'^^^ 



and AV:^^Zt 
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iVesperini. Zepf fc Kundu |l2003j) do not shape the Som- 
brero GCLF in an observable manner. 

3.3. Blue Tilt 

A visual inspection of the color-magnitude diagram in 
Figure [Till suggests a correlation exists between color and 
magnitude for metal poor GCs, while no such relation 
is apparent for the metal-rich GCs. Quantifying this 
"blue tilt" requires careful consideration of the biases 
due to parameter limits, photometric errors, and con- 
taminants. For instance, the introduction of a subpopu- 
lation boundary at a single color (i.e. blue GCs strictly 
have B— R < 1.3) will bias an examination of the indi- 
vidual subpopulation color-magnitude distributions. For 
this reason, the Nmix routine (introduced in W6.ll is used 
to quantify the subpopulation color peaks at five mag- 
nitude intervals. Two lines were then fitted to the blue 
and red Nmix peaks with a least-squares minimization 
routine that was weighted by the Nmix peak standard 
deviation in the mean. Weighting the color peaks by 
this value effectively propagates photometric error infor- 
mation, since less weight will be placed on the higher 
color dispersion characteristic of fainter GCs. 

The magnitude intervals were chosen to cover the en- 
tire luminosity range while still maintaining a numerical 
significance at each interval. Since the faint end of the 
magnitude distribution is more affected by photometric 
uncertainty, the faintest magnitude interval is excluded 
from the final linear fit but is shown in Figurc lTUI Indeed, 
whether the faintest bins follows the trends established 
by the brighter points is indeterminate. Table 01 sum- 
marizes the slopes from linear fits to all band and color 
combinations. 

The formal significance of the linear slopes from the 
B— R Nmix peak fits are consistent with the visual im- 
pression that a more significant slope is present for the 
blue GCs, 5(B-R)/(5R -0.035 ± 0.005, compared to 
the red GCs: S{B-R)/5R = -0.015 ± 0.007 (see the fits 
in FigurelTUIl. From Tabled it is evident that for a given 
color, the slope (and the formal significance) changes de- 
pending on which band is used on the abscissa. When 
the midpoint between the red subpopulation (5(B— R) / (5B 
and 5(B-R)/(5R slopes is considered, -0.012 ± 0.007, it 
provides no strong evidence for an analogous red GC 
trend. 

Two recent ACS studies of GC systems discovered 
similar blue t ilts in the metal-poor subpopulation. 
IStrader et all l|200(iD showed that the massive Virgo 
cluster galaxy, M87 (and perhaps NGC 4649, but not 
NGC 4472), has a metal-poor color-magnitude trend in 
the GCs brighter than GCLF TOM. The Sombrero trend 
extends to GC at least V ^ 22.7 (My ^ —7.1), consistent 
with the Sombrero GCLF TO M and the mag nitude range 
found in lStrader et al. Il|2006|) . iHarris et alTl)2006D found 
similar trends in eight brightest cluster/group galaxies. 
However, they suggested that the blue tilt is found only 
among the brightest of the blue GCs and gradually be- 
gins to disappear at a transitional magnitude approxi- 
mately 1.5 mag brigh ter than the GCLF TOM. Neither 
IStrader et all l)200m nor 'H arris et all lj2006|) found any 
evidence for an color-magnitude trend among the metal- 
rich GCs. The discovery of a blue tilt in a field spiral like 
Sombrero shows this phenomena is not restricted to the 
GC systems of massive ellipticals in high-density galaxy 



environments. Blue tilts have since been discovered in 
o ther spiral ga l axies ( Forde et al. in prep.). 

IHarris et al.l l)2006|) converted their color-magnitude 
trend to a metal-poor GC meta llicity-mass gradient: 
Z oc M°-55. IStrader et aD ^200(f] found similar values 
of Z cx in M87 and Z oc M"-^^ in NGC 4649. 

The corresponding gradient from the Sombrero data, 
Z oc M"-^^, is from the midpoint of the metallicity- 
mass proportionalities from the 5B—K/SB and 5B—K/SR 
color-magnitude slopes. Systematics from band and 
color-to-metallicity transformation differences cannot be 
ruled out as the cause for the smaller Sombrero gradi- 
ent. However, a similar color-to-m etallicity transforma- 
tion technique to IHarris 'eran i|20061 was used here, sup- 
porting a real difference in the trends. 

Two interpretations of the GC metallicity-mass trend 
were proposed in the earlier studies: "pre-enrichment" , 
where larger proto-GC clouds contain higher metal abun- 
dances and "self-enrichment" where GCs were able to 
retain and then incorporate some metals created within 
the GC. Harris et al. (2006) supported a pre-enrichment 
origin after noting that some numerical simulations pre- 
dict that only the most massive dE galaxies (> lO^M©) 
can retain metals, while anythi ng less massive will 
tend to eject mctal-rich gases fe.g;.. lMac Low fc Ferraral 
^ 99: Fragile Murr av & Lin 20Q3)- It was suggested by 
Harris et alJ ijiOOfil) that dEs in this mass regime are ca- 
pable of building GCs with masses well above IO^Mq, 
a value which corresponds to slightly more than half a 
magnitude brighter than their transitional magnitude of 
Mi = -10.0 (they assum e M/L = 3). The Sombrero 
and IStrader et al~1 l)2006(l results support a fainter tran- 
sitional magnitude (~ 1.5 mag fainter) indicating the 
trend is apparent among objects less mass ive than the 
mass limit for GC pre-enrichment quoted in IHarris et alJ 
(2006). 

In terms of absolute metallicity, metal-rich GCs have 
on average 10 times higher metallicity than metal- 
poor GCs, thus any given metallicity-magnitude gradi- 
ent found among metal-poor GCs would correspond to 
a 10 times smaller slope for the metal-rich GCs, which 
may have disappeared in the noise. If the metal-rich 
GCs actually do not follow an analogous trend, differ- 
ent formation environments (e.g. Strader et al. 2006 and 
Harris et al. 2006) or subpopulation characteristics (e.g. 
differing radial distributions) are both conceivable ideas 
that would require careful theoretical modelling to test. 
Indeed, without predictions from modelling of individual 
GCs formation and evolution, attempts to explain physi- 
cally the processes responsible for this new discovery will 
largely remain speculative. 

3.4. Radial Distribution 

The Sombrero metal-rich GCs are more centrally- 
concentrated than their metal-poor counterparts, as 
shown in Figure 1141 From Figure ^ it is visu- 
ally apparent the metal-rich GCs are associated with 
the Sombrer o bulge rather than the disk, confirming 
the work of iMoretti et aO ll2003D and suggestions by 
iForbes. Brodie fc LarsenI l)200H) . The surface densities 
of the entire GC system and the two subpopulations show 
slopes which are similar to the galaxy light profile (Fig- 
ure^)) with neither subpopulation displaying a particu- 
larly better match. 
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Fig. 14. — The GC surface densities of the entire population 
(crosses) and the two subpopulations (square and triangle for blue 
and red GCs, respectively). The solid, dotted, and short-dashed 
lines are de Vaucouleurs profile fits to the entire, blue and red sur- 
face densities, respectively. The lower panel is a log representation 
while the top panel shows a linear view. In the upper panel the 
dash-dotted line and circle symbols are from ground-based analysis 
of the Sombrero GC surface density (Rhode & Zepf 2004) and the 
long- dashed line is from HST WFPC2 (Larscn, Forbes & Brodi<3 
120011) analysis. Error bars are from Poisson statistics. For clarity, 
the online version of this Figure is presented in color. 
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Fig. 15. — Surface densities of the blue subpopulation (top), red 
(middle) and the entire GC system (bottom). These are compared 
with the arbitrarily normalized R-band galaxy luminosity profile. 
The GC surface densities all resemble the galaxy profile. 



Two functions traditionally used to model 
GC surface densities {(Jgc) are a power-law, 
log(CTGc)= apo«)-log(Rgc)+/3po«), and de Vaucouleurs 

law, log(crGc)= adev • Rgi^ + PdeV- A shallower 
power-law exponent, otpow = 

-1.14 ± 0.13, is found 
compared to the surface density exponent derived 
recently from ground- based, wide-field imaging: 
apow ^ -1-85 ± 0.07 lIRhode fc Zeofl l200l . The 



de Vaucouleurs fit also yields a shallower curve of 
c^dev — —1.77 ±0.17 compared to ground-bas ed and 
HST WFPC2 ijLarsen. Forbes fc Brodiel l200l|) anal- 
ysis, which found slopes of adeV = —2.11 ±0.08 and 
adeV = —2.00 ± 0.18, respectively. The masked Som- 
brero disk region was excluded from the surface density 
calculation. 

The current ACS and past WFPC2 de Vaucouleurs 
surface density fits to the entire GC population are com- 
pared in Figure PPil Given that the WFPC2 fit was de- 
rived from only three pointings (thus inherently uncer- 
tain due to both the small areal coverage and GC sample 
size), it is likely consistent with the ACS results within 
the uncerta inty due to sampling statistics. While the 
fit from the iRhode fc Zepf I 12004*) is clearly offset from 
the ACS curve, the surface density points from this past 
study are consistent with the ACS fit. This supports 
the idea that a shallower surface density is present in 
the central galaxy regio ns compared to regions outside 
of the ACS field of view. Bridges fc Ha^lTool found 
that Sombrero profile could be best represented by a 
composite of two power laws transitioning at 2.5' with 
a shallower profile in the inner radial region. This is also 
consistent with observations of other GC systems that 
exhibit flattened cores in their GC radial distributions 
(e.g. Dirsch ct al. 200$. 

Adopting an outer radius of 19' (where ace ^ 
arcmin"^; Rhode & Zepf 2004) and integrating the de 
Vaucouleurs surface density profile, a total GC popu- 
lation of 1870 ± 850, a blue total of 1110 ± 850 and a 
red total of 820 ± 200 were estimated (population errors 
are la values from Monte-Carlo simulations of the inte- 
gration of the surface density profile). The total pop- 
ulation deriv ed here is consi s tent w ith the 1900 GCs 
estimated bv IRhode fc Zeofl l)2004j) . However, the de 
Vaucouleurs profile continues to increase significantly as 
Rgc ^ and may result in an overestimation of the 
total. Indeed, the de Vaucouleurs profile of the entire 
GC population overestimates the number of GCs found 
in the innermost radial intervals {]<f^^'^'^*'^'^ ~ j^3q jj^ 
0.3' <R(,c < 1-0') compared to the number of observed 
GCs {Wd^^''''^'^ « 90) by « 30%. Fitting a "cored" pro- 
file of the form ace = a • (1 + {-Cgc/l'YY (from the 
empirical surface density fit in Dirsch et al. 2003; see Fig- 
ure^)) to the surface densities of entire population, and 
the blue and red subpopulations separately, yields the re- 
duced estimates of 1540 ± 640, 940 ± 620, and 670 ± 270, 
respectively (errors computed analytically). 

Forbes, Brodie fc Larscn (2001) showed that the num- 
ber of red GCs normalized to the bulge luminosity, the 
bulge specific frequency (bulge Sat), has a similar value in 
the Milky Way, M31 and the Sombrero, perhaps indicat- 
ing a common bulge formation history. To compute the 
bulge Sat a fraction of the total galaxy luminosity must 
be attributed to the bulge component. The Sombrero 
reddening-corrected total By magnitude is 8.39 ± 0.06 
and the corrected color is B-V= 0.87 ± 0.01 (RC3), 
which corresponds to My = —22.25 ± 0.06 for a dis- 
tance modulus of 29.77± 0.03. Adopting a bulge-to-total 
luminosity ratio of 0.73 ijBaggett. Baggett fc Anderson I 
1998), resulted in a bulge apparent magnitude of My = 
-21.91 ±0.06. 

A bulge Sat = 1.2 ± 0.8 was calculated (errors from 
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Fig. 16. — The GC surface d ensit ies. Same as Figure IT^ except 
a "cored" profile (presented in H'iAt is fit to the surface densities. 
This profile provides a more accurate population estimate because, 
unlike the de Vaucouleurs profile, it does not continue to increase 
significantly as Rgc — * 0. 



the GCs totals) from the extrapolated "cored" red sur- 
face density fit (N^ed — 670 ± 270). This value is 
larger than th e bulge Sn estimated for M31, the Milky 
Way l|Forbes. B rodic fc Larson 2001), and five spiral 
gal axy GC systems analyze d in a WFPC2 compilation 
by iGoudfrooii et al. I l)2003|) . In these previous stud- 
ies a higher degree of uniformity was found if instead 
a bulge Sn was calculated from the number of red 
GCs within two effective radii (R1/2) of the galaxy 
bulge. Adopting a Sombrero bulge effect ive radius of 
0.89' l|Baggett. Basfgett fc Anderson II1998D . 160 ±30 red 
GCs within 2R1/2 are estimated, which corresponds to 
a bulge Sn = 0.3 ± 0.05. This version of the Sombrero 
bulge Sm IS consistent with M31 (bulge Sn = 0.4 ± 0.2; 
iForbes. Brodie fc L arsen 2001) and is larger than all of 
the spirals studied in .Goudfrooii et al. 1 1)2003|) . 

The M31 bulge Sat is likely more uncertain than the 
formal errors quoted above and an overestimate since the 
number of metal-rich G Cs is derived mo stly from photo- 
graphic plate analysis ijFusi Pecci et alJfT993|) . which is 
poorly suited to address the significant amount of con- 
tamination expected because of the face-on orientation 
of the M31 disk. Thus the Sombrero bulge Sn value is 
likely larger than all the spiral GC systems on which this 
analysis has been performed. If spirals generally have a 
constant bulge Sat, it would support the idea that most 
metal-rich GCs formed in conjunction with their host 
bulge and the Sombrero analysis would seem to suggest 
that something characteristic of more massive galaxies 
(e.g. mass-to-light ratios, formation timescales) leads to 
enhanced metal-rich GC formation and/or decreased de- 
struction rates. 

When GC Nmix peak colors derived from GCs brighter 
than V = 22.0 are plotted against Rgc (Figure no 
significant trend is observed in the red (—0.000 ± 0.015 
mag arcmin"^) and blue subpopulation (—0.014 ± 0.012 
mag arcmin^^). Figure [T7I also shows that the galaxy 




1 - 



1 2 3 

Rg^(arcmin) 

Fig. 17. — GC color versus projected galactocentric radius. 
Crosses are those GCs with V < 22 magnitudes. No statistically 
significant gradient is found in either the red (dashed line and tri- 
angles) or blue subpopulations (dotted line and squares). Error 
bars are the standard deviation in the means of the peak color es- 
timates. The mean B— R values for the total GC population (solid 
line and stars) illustrate the changing ratio of blue to red GCs with 
radius. The thick solid line shows a fit to the ACS mosaic galaxy 
color profile. 



color profile extracted from the ACS mosaic is redder 
than the metal-rich subpopulation but begins to show 
GC-like colors at ^ 2'. 

3.5. Sizes 

Figure IHl hints at a size-color trend for the GCs. His- 
tograms of the half-light radii of the blue and red GC 
subpopulations are presented in Figure 15. The blue 
and red subpopulations have median R/i/ of 2.12 pc and 
1.87 pc, respectively; thus the reds are systematically 
~ 13% smaller than the blue GCs (mean sizes 2.40 pc 
and 2.10 pc; reds ~ 14% smaller). Using the same 
size determination algorithm (ISHAPE) and a Sombrero 
WFPC2 F547M-band image centered on the galaxy cen- 
ter, iLarsen et alTl |2001) found mean half-light radii of 
2.16 pc and 1.67 pc (updated with the current distance 
modulus) for the blue and red subpopulations, respec- 
tively, corresponding to a larger size difference of about 
30%. If the ACS size measurements are restricted to the 
GCs found in the same region as the central WFPC2 
image (see Figure ^ , a blue-red mean size difference of 
30% is found, with mean R^,; of 2.25 pc and 1.73 pc for 
the blue and red GCs, respectively. This size difference 
is consistent with the WFPC2 data, suggesting that the 
smaller average size difference for the full ACS mosaic is 
a real effect, due to the larger area covered by the ACS 
data. 

Two main interpretations of the size difference between 
blue a nd red GCs have been proposed. Larsen & Brodie] 
(|2003) suggested typical observations of the central com- 
ponent of a GC system could lead to a subpopulation 
size difference from the projection of the different GC 
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Rh,(pc) 

Fig. 18. — GC half-light radius histograms for blue (top panel) 
and red (bottom panel) GC subpopulations. The red GCs appear 
13% smaller on average than the blue GCs with median R^j values 
of 2.12 and 1.87 pc for the blue and red subpopulations, respec- 
tively. 
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Fig. 19. — GC sizes versus galactocentric radius normalized 
by the Sombrero bulge effective radius. Left, middle, and right 
sections contain the entire population, red subpopulation, and blue 
subpopulation, respectively. Solid lines are power-law fits to the 
entire radial coverage while the dashed lines are fits only to the 
inner GCs (R9C/R.1/2 < 3). The difference between the fit to the 
inner and entire blue GC sample suggests projection effects are 
more apparent in the central regions for this subpopulation. 



subpopulation radial distributions (with the red GCs be- 
ing more centrally-concentrated) and an underlying GC 
size dependence on the three-dimensional galactocentric 
radius (Ran) of all GCs. 

Alternatively, in iJordan I 1)20041) it was shown that a 
correlation of GC half-light radii and metallicity can arise 
from mass segregation combined with a metallicity de- 
pendent main-sequence lifetime. In the same study, it 
was shown that a form of this model fit observations 
of M87 GCs. It is import ant to n o te th at the size- 
color trends predicted by the IJordan I l)2004j) models are 
strongly dependant on the relative ages of the two sub- 
populations. To reproduce observations IJordan I 1)20041) 
assumed the subpopulations are coeval. While GC sys- 
tem formation scenarios generally place metal-rich GC 
formation later than the metal-poor GCs and the ages of 
Milky Way GC may show an age-metallicity relationship 
(e.g. (Salaris fc Weiss 200^), a precise spectroscopic con- 
strai nt on extragalactic GC ages is not currently achiev- 
a ble llStrader et al. Il2005a|) . 

IJordan et'^Lri(2005)l addressed both of these theories in 
their analysis of a large sample of GC systems from early- 
type galaxies in the ACS Virgo Cluster Survey (VCS; 
iCote et al. 11200^ . They concluded the intrinsic metal- 
licity difference between blue and red GCs provides a 
better explanation for the observed GC subpopulation 
size differences than projection effects largely because the 
size-Rgc trends they observed were suggested to be too 
shallow to cause the observed mean size difference. How- 
ever, the size-Rgc trends they examined were restricted to 
the metal-poor subpopulation, which is more affected by 
projection effects compared to metal-rich GCs as demon- 
strated below. 

Irrespective of the previous results, IJordan I l)2004l) 
noted the only way to unequivocally distinguish between 



the two theories is to examine GCs at large projected 
radii, since the GCs furthest from the projected galaxy 
center will be more likely to share a similar R315 value, 
thus enabling a comparison between the subpopulations 
without the influence of an underlying size-Rsu trend. 

The ACS mosaic coverage extends to w 6 times 
the Sombrero bulge effective radiu s (R1/2 = 0.89'; 
iBaggett. Baggett fc Anderson 1 1199^. T his is roughly 
twice the radial limit reached bv I Jordan e t al. ( 2005) in 
their most massive galaxies, whose populous GC systems 
constituted the bulk of their size sample. FigurelT^shows 
the clear GC size-(Rgc/R3D) trends in the Sombrero sub- 
populations and the entire population. When a power- 
law relation of the form log(R;i/)= a + 6-log(Rgc/Ri/2) 
(as in IJordan et al.l 12005)) is fit to the data, the fol- 
lowing slopes are derived: huue — 0.16 ± 0.04, bred — 
0.32 ± 0.05, and = 0.24 ± 0.03. For comparison, 
IJordan et alJ l)2005)) found for the metal-poor GCs an 
average weighed slope of = 0.07 ± 0.01 from their 

sample of Virgo GC systems. Fitting curves to the Som- 
brero size-Rgc information restricted to the Virgo radial 
limit (Rgc/Ri/2 < 3) yields shallower exponents com- 
pared to the unrestricted ACS values and a blue subpop- 
ulatio n value which is consistent with the IJordan et alJ 
1)2005)) value: buue = 0.09 ±0.07, bred = 0.27 ±0.07, and 
6an = 0.19 ±0.05. 

Figure[2n|effectively contrasts both the GC subpopula- 
tion size theories and clearly illustrates that there is little 
subpopulation size difference at larger Rgc. This result 
was alluded to at the beginning of this section when the 
ACS sizes were compared to the results from the Som- 
brero WFPC2 image, which was centered on the galaxy 
(see Figure ^1 . The nearly flat blue GC size-Rgc gradi- 
ent in the inner region (< 3 Rgc/Ri/2; Figure lT9)l could 
be caused by the introduction of scatter from the wider 



14 



Spitler et al. 




1 1.2 1.4 1.6 1 1.2 1.4 1.6 10 a 6 10 a 

B-R M^, M„ 



Fig. 20. — Size-color plots from two equally-sized galactocentric 
radius intervals divided at R^c = 2.4'. The solid lines represent 
moving biweight location estimators from 50 data values with la 
boundaries shown. In contrast to the inner radial interval (left 
panel), the biweight values from the outer interval (right panel) 
are consistent with no size-color trend. The percentage difference 
between the subpopulation sizes is 17 ±5% and 5 ±6% for the inner 
and outer intervals, respectively. These observations are most likely 
a result of projection effects and the different radial distribution of 
the two subpopulations. 



Fig. 21. — Globular cluster subpopulation half-light radii versus 
magnitude. The line is a moving biweight location estimator (width 
of 50 data points) with Icr boundaries. The metal-poor GCs show 
no clear size-magnitude trend among the fainter GCs until < 
—9.0 mag where a strong upturn is observed. The metal-rich GCs 
are consistent with this upturn and also show increasing sizes with 
luminosity along the entire magnitude range. 



range of blue GC R^d values (i.e., projection effects). 
These observations suggest projection effects are respon- 
sible for the appearance of a Sombrero subpopulation size 
difference. 

To conclude, the observed Sombrero s ubpopulation 
size di fference is consistent with the data in lJordan et al.l 
(j2(},Q5|) , although the extended radial coverage leads to 
a different interpretation. Projection effects are most 
likely responsible for the observed average size difference 
between blue and red GCs in the Sombrero and that, 
in general, projection effects should always be accounted 
for when comparing GC sizes from different GC systems. 
The Sombrero size analysis reaffirms that GC size in- 
formation at large Rgc is essential for disentangling the 
two interpretations of GC subpopulation size difference. 
Also, it is important to examine the metal-rich subpop- 
ulation for size-Rgc trends since this subpopulation is 
less susceptible to projection effects from its centrally- 
concentrated distribution. 

3.6. Size- Magnitude Trend 

Figure |2] shows the size- magnitude diagrams of both 
GC subpopulations. Solid lines show moving biweight 
location estimator values from the ROSTAT software 
routine ffieers et al.lll990|) with a window width of 50 
data points. Among the metal-poor GCs, no clear size- 
magnitude trend is detected from the faintest GC mag- 
nitudes to My — —9.0 where a strong size-magnitude 
"upturn" is observed such that the brighter GCs show 
increasingly larger sizes on average. The metal-rich GCs 
exhibit a shallower size-magnitude trend over the entire 
magnitude range where the average sizes again become 
larger as the GC luminosity increases. 

In i)2.1.2l it was shown that the la uncertainties from 



a linear fit to the size-size plots derived from each of the 
three image bands corresponds to an accuracy level of 
^ 0.2 pc (FigureEJ, which is approximately half the dy- 
namical range of the size-magnitude trends. When the 
same analysis is performed on GCs of three different mag- 
nitude intervals, it is estimated that the sizes are accu- 
rate at the 0.15, 0.17, and 0.27 pc levels for the brightest 
(My < -8.3), intermediate (-8.3 < My < -7.4) and 
faintest (My > —7.4) thirds of luminosity distribution, 
respectively. 

To determine the whether the shallow trend is signifi- 
cant, the two subpopulations and the entire sample were 
bootstrapped and divided into magnitude intervals of ap- 
proximately 50 data points. The biweight of each interval 
was calculated and a least-squares linear fit performed 
on the resultant biweight values for all but the brightest 
magnitude interval to prevent the size-magnitude upturn 
from dominating the fit. Figure [521 contains histograms 
from the output of this procedure and suggests a sig- 
nificant trend is present among the metal-rich GCs and 
marginally so in the metal-poor subpopulation. 

3.6.1. Size- Magnitude Upturn and More Massive 
Stellar Systems 

To date, no positive size-luminosity trend has 
been reported in aii y extra galactic GC s ystems (e.g. 
IKundu fc WhitT^IT;;^ iOOl'; ■T.aTseu et a,l. I ^Q^, al- 
though there exists no comparably large dataset of GC 
sizes from such high signal-to-noise images as is pre- 
sented here. There are, however, a growing number 
of studies in which very luminous and extended GCs 
are being found. For example in NGC 5128 four- 
teen massive (My > -9.6; M > 10*^ M^) GCs were 
shown to ha ve unusually large sizes (Martini & Ho 200^ 
iHarris et al. 2002) . Martini fc Ho . 1,2004,) hypothesized 
these to be accreted nucleated dwarf elliptical (dE,N) 
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Fig. 22. — Histogram of slopes from a bootstrap statistical anal- 
ysis of the size-magnitude trends for the blue, red, and entire GC 
sample in the bottom, middle and top panels, respectively. The 
high-luminosity size-magnitude upturn was intentionally avoided 
in this analysis, as described in the text. This analysis suggests a 
significant size-magnitude trend exists among the metal-rich GCs 
and marginally so among the metal-poor GCs. 



a. 




Fig. 23. — Globular cluster half-light radii versus magni- 
tude. The line is a moving biweight location estimator (width 
of 50 data points) with Icr boundaries indicated by the lines. 
The dashed line represents the division between Milky Way GCs 
and the susp ected Local Group accreted dE,N nuclei candidates 
IMackev Ik van den Bergh ,,2005.') . 



nucl ei that were stripp ed of their halos. Using ACS 
data IHarris et all f2006) found 20 to 30 bright and ex- 
tended objects in NGC 1407. Also certain Galactic GCs 
brighter than My = -9.4 (uj Gen, M54, and NGC 



2419) are all larger than the median Milky Way GC 
size and are suggested to be accreted dE N nuclei (see 
refs. in IMackev fc van den Bergh 1 12005|) . Figure OS 
shows that the boundary between those suspected Local 
Group accreted dE,N nuclei and typical Milky Way GCs 
iMackev fc van den Berghl [2005 ) approximately corre- 
sponds to where the size-magnitude upturn appears in 
the Sombrero. 

These observations support the idea that a group of 
particularly massive and extended GCs may exist in 
other galaxies, indicating the Sombrero high-luminosity 
GC size-magnitude upturn may not be unique. It is also 
apparent that utilizing an accreted dE,N nuclei scenario 
to explain the presence of massive and extended GC- 
like objects is common practise. This idea is supported 
by a resemblance between GCs (metal-poor GCs in par- 
ticular) and nuclei still embedded within their dE,N 
galaxy. Similarities between massive, met al-poor GCs 
and dE nuclei include: broadband colors lILotz et al. I 
I2OOI IDe Propris eTaTI l2005t IHarris et al.ll2006|) . the 
peak in the dE n uclei luminosity fu nction {My ~ —9.5; 
from table 2 of iLotz et al. 1 1200^ corresponds to the 
most luminous GCs (and where the Sombrero GC size- 
magnitude upturn appears), and dE nuclei and metal- 
poor GCs show a nearly identical color-magnitude trend 
(Harri s ct al 2006) . 

iiCote et al.l l)200fil) studied the nuclei of 51 early- type 
galaxies from the ACS VCS and showed that the sizes 
of the nuclei increase with their luminosity. If the nuclei 
size-magnitude trend is partially preserved during the ac- 
cretion of a dE,N onto the Sombrero, a size-magnitude 
trend might begin to appear as the number of accreted 
nuclei becomes comparable to the high- luminosity indige- 
nous metal-poor GCs. 

Another possible explanation for the size-magnitude 
upturn is that it reflects an intrinsic property of GC for- 
mation. This idea is motivated by an analogous size- 
mass upturn found in t he young counterparts t o GCs, 
young massive clusters ()Kissler-Patig et al.ll2006l) . Sup- 
porting this is a striking absence in Figure 1231 of Som- 
brero GCs with My < —10 that are smaller than the 
median size of the sample (Rm ~ 2 pc). This is unex- 
pected given the intrinsic scatter in sizes is quite large 
(Sombrero size standard deviation is 1 pc) and if the 
upturn was caused entirely by the presence of accreted 
dE nuclei, high-luminosity indigenous GCs with small 
sizes should still be observed. Thus an absence of small 
bright objects is consistent with an intrinsic origin for 
the size-magnitude upturn. 

Similarly determining if a deficiency of bright GCs with 
small sizes exists in other systems (and therefore a size- 
magnitude upturn) is a simple exercise to test whether 
the size-magnitude uptu rn is a common feat ure. An ex- 
amination of figure 6 in iMieske et alJ l|200(iD containing 
the large GC size-magnitude dataset from the ACS For- 
nax Cluster Survey (Jordan et al., in prep.) suggests a 
such deficiency exists in other GC systems. 

[Mieske et ah, (^2QQ6) proposed that all compact stellar 
systems (in particular Ultra Compact Dwarfs, UCDs) 
with My < —11.0 show increasing sizes with increas- 
ing luminosity, while less luminous stellar systems have 
magnitude-independent sizes and therefore are likely to 
be typical GCs. The figure that partially motivated this 
idea is the same which here is suggested to demonstrate a 
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deficiency of iuminous objects witti smali si zes. TIius tlie 
size-m agnitude upturn shown in figure 6 of iMieske et al.l 
()2006j) may actually begin at My ^ —10.0 where the de- 
ficiency noted here becomes apparent. A fainter absolute 
magnitude marking the emergence of the Sombrero up- 
turn, Mv ^ —9.0, may indicate the upturn appears in 
even less massive compact ste llar systems than demon- 
strated in IMieske et al.l l)2006() . If the Sombrero upturn 
is indeed an intrinsic feature of its massive GCs, it sup- 
ports a continuous trend of magnitude-dependent sizes 
from the brightest GCs to more massive stellar systems 
such as dE nuclei and UCDs. 

3.6.2. Metal- Rich GC Size- Magnitude Trend 

Validating the shallower metal-rich size-magnitude 
trend is difficult because completeness effects (where 
faint, extended GCs are missed) conceivably could have 
caused the trend. A quantitative investigation into these 
effects and a careful consideration o f the numerous size- 
related trends (as demonstrated in I.Tordan et alJl200fi^ 
is beyond the scope of the present work. It is also noted 
that a combination of the extended radial distribution of 
the metal-poor subpopulation and an underlying size- 
Rgc trend (i.e., projection effects, see i l3.5|l plausibly 
could introduce enough scatter to hide an analogous size- 
magnitude trend in the metal-poor subpopulation. 

Even so, similar size-magnitude trends have been found 
among young star clusters of the merger galaxy NGC 
3256 (Zepf et al. 1999) and 18 nearby spiral galaxies 
ijLarsenI 12004(1 . The former study found that the clus- 
ters follow a size-luminosity proportionality of Hhi oc 
^~0.07^ compatible with the results of the later spi- 
ral study. The corresponding Sombrero metal-rich GC 
fit, Kfii oc ^o.o7±o.03^ jg identical to the past results. 
A fit to biweight location estimates of five Sombrero 
metal-rich GC magnitude intervals confirm this analy- 
sis: Rhi oc l0 09±o.oi_ ^Yjjiig jg curious that the results 
are consistent across the different studies, it is still in- 
determinate whether the consistency is a reflection of a 
real trend or the result of a similar type of completeness 
effect. 

4. SUMMARY AND CONCLUSIONS 

A six-pointing mosaic from the HST Advanced Cam- 
era for Surveys was used to study the globular clus- 
ter system of the Sombrero galaxy. An unprecedent- 
edly deep 95% of GCLF and ^ 2 mags fainter 
than the turnover magnitude), virtually contamination- 
free sample, allowed an accurate GC luminosity function 
(TOM) turnover magnitude to be determined: iHy^^^ — 
22. 17 ±0.06 or M^^*^ = -7.60 ± 0.06, assuming m-M= 
29.77±0.03. The GCLF was found to be well fit by a tg to 
2 magnitudes below the TOM fFigure [TT|l . The disper- 
sion of the t5 fit is CTj^^-^ = 1.08±0.06. While the TOM is 
slightly fainter t han a universal TO M, M^°*^ = -7.35 
l|Richtler II200I I Jensen et al. \\200^ . it is likely consis- 
tent within the typical observational errors. No variation 
with Rgc is found in the Sombrero GCLF TOM, nor is a 
significant difference between the subpopulation TOMs 
detected. 

The GC color distribution is clearly bi-modal (Fig- 
ures 0] and with blue and red peaks at B— R = 1.15 
and 1.44 or [Fe/H] = -1.38 and -0.49, respectively (see 



Table 121). Subpopulation peak colors are consistent with 
the expected values for a galaxy with a luminosity like 
the Sombrero, according to the relation between the GC 
subpopulation color peaks and host galaxy luminosity 
l)Strader et al.l 120041. It was found that the color dis- 
persion in the metal-poor subpopulation is smaller com- 
pared to the metal-rich subpopulation. This is in agree- 
ment with the results of Harris et al. (2006), Peng et al. 
(2006), and Strader et al. (2006). 

An estimated 53% of this GC sample belongs to the 
blue subpopulation. This proportion contrasts to the re- 
sults from ground-based observations where a larger pro- 
portion of 59% was found from data that covers the entire 
spati al distribution of the GC system (Rhode fc Zepf I 
|2004() . Thus the difference is likely caused by the decreas- 
ing ratio of red to blue GCs with galactocentric radius 
fFigurc [T7|l . No color-Rgc trend is observed within either 
subpopulation. The galaxy color profile tends to be red- 
der than metal-rich GCs and only begins to overlap this 
subpopulation at Rgc > 2'. 

The best-fitting de Vaucou leurs radial pro fi le of the 
Sombrero GC system from iRhode fc Zepf I l|2004f) is 
steeper than the one found here. However, their sur- 
face density values from within the ACS mosaic region 
are consistent with the ACS surface densities (FigureEJ. 
This observation is consistent with there being a core of 
constant number density in the inner GC system, con- 
firming observations m ade in a n earlier Sombrero GC 
system study (Bridges fc Hanesl [l992'l. 

The number of red GCs normalized to the bulge lumi- 
nosity, the bulge specific frequency, is found to be larger 
than all spiral galaxies for which this value has been cal- 
culated. This suggests the rate of metal-rich GC for- 
mation was enhanced and/or GC destruction was less 
effective in the Sombrero galaxy relative to less massive 
spirals. 

A color-magnitude trend, a so-called "blue tilt" , was 
discovered in the metal-poor GC subpopulation that is 
likely a metallicity-mass trend (Figure ^J. The metal- 
rich subpopulation shows no strong evidence for an anal- 
ogous trend. The Sombrero provides the first example 
of this trend in a spiral galaxy and a galaxy found in a 
low-density galaxy environment. The blue tilt extends 
from the brightest GCs to at le ast the GCLF TOM , 
in agreement with the findings of iStrader et al. I l)2006j) . 
A shallower trend, in terms of metallicity-mass propor- 
tionality, is found here compared to examples in mas- 
sive ellipticals found i n high-d ensitv galaxy environments 
l|Strader et al. II2006I: iHarris e t al. 2006). This difference 
and the absence of a blue tilt in NGC 4472 likely indi- 
cates the gradient of the color-magnitude trend depends 
on the specific properties of the host galaxy. 

It is shown that GC sizes correlate with Rgc (Fig- 
ure 119(1 . From the extended radial coverage provided 
by the Sombrero ACS mosaic, the metal-rich GCs are 
found to have a steeper size-R^c gradient compared to the 
metal-poor subpopulation, thus causing the mean sub- 
population sizes to become identical at larger Rgc (Fig- 
ure 120)). These observations are consistent with projec- 
tion effects being largely respo nsible for the observ ed GC 
subpopulation size difference (|Larsen fc Brodiel 12003). 
as opposed to an intrinsic feature of their respective 
metallicities (^Jordan 2 0041 . 

A size-magnitude upturn was discovered among the 
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brightest GCs (My < —9.0), where the sizes increase 
as a function of increasing luminosity (Figure [2J- The 
accretion of dwarf elliptical nuclei into the Sombrero GC 
system provides a plausible scenario for this trend. How- 
ever, a n analo gous trend among young massive clusters 
iKissler-Patig et al.l l)200(i() and a noticeable deficiency of 
luminous GCs (My < —10.0) with smaller than average 
sizes (Figure 1231 also see figure 6 in Mieske et al. 2006), 
supports an intrinsic origin for the size-magnitude up- 
turn. An intrinsic size-magnitude upturn suggests the 
most massive GCs are at the low-mass end of a contin- 
uous trend that includes more massive systems such as 
dwarf nuclei and UCDs. 

A shallower positive size-magnitude trend among 
metal-rich GCs is found at all magnitudes. The 
luminosity-size proportionality, Khi oc ]^^-07±o.03^ jg 
sistent with the same fits to young star cluste rs of 
the merger galaxy NGC 3256 llZepf et alJ 11999(1 and 
18 nearby spiral galaxies l|Larsenl 12004(1 . It is unclear 
whether the consistency is a reflection of a real trend or 
the result of a completeness effect. 
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TABLE 1 

Aperture and Reddening Corrections 



Filter F435W (B) F555W (V) F625W (R) 



AC(5-10 pixels) 


0.187 


0.191 


0.198 


AC(lO-infimty) 


0.107 


0.092 


0.088 


Total AC 


0.294 


0.283 


0.286 


Mean Extinction Correction (Ax) 








from the six pointings 


0.210 


0.163 


0.136 



TABLE 2 

Color/Metalicity Distributions and Subpopulation Proportions 







Blue 


Red 


Blue [Fc/H] 


Red [Fc/H] 




^red 


JFe/H] 
blue 


JFe/H] 
red 


Blue % 


B- 


~R 


1.15 


1.44 


-1.32 


-0.46 


0.08 


0.10 


0.23 


0.29 


53% 


B 


-V 


0.67 


0.90 


-1.41 


-0.50 


0.05 


0.07 


0.24 


0.31 


52% 


V- 


-R 


0.46 


0.55 


-1.24 


-0.52 


0.02 


0.03 


0.22 


0.32 


52% 



TABLE 3 

Sombrero Distance Moduli Estimates 



Technique (source) 


m-M 


Mpc 


Notes 


Surface Brightness 


29.79 ± 0.18 


9.1 ±0.8 


Re-calibrated by 


Fluctuations (SBF) 






—0.16 mag. 


(Tonrv et al. 2001) 






C Jensen et al. 2003") 


Fabcr- Jackson 


30.03 


10.1 


Bt = 9.41 ±0.02 


relation 






cr = 237 ±6 


(HyperLeda) 






m-M=BT + 6.21ogo- + 5.9 


Virgo infall 


30.85 


14.8 


Close proximity makes 


velocity 






velocity uncertain 


(HyperLeda) 






Ho = 72km/s/Mpc 


Planetary Nebulae 


29 74+U-U4 


8.9 ±0.6 




LF (PNLF) 






(Ford et al. 1996) 








Adopted distance 


29.77 ±0.03 


9.0 ±0.1 


Weighted average of 


to Sombrero 






SBF & PNLF 




TABLE 4 




GC Color-Magnitude and Metallicity-Mass 




Gradients 




color /band 


Blue Color 


Rod Color Z oc M°' 


(B-R)/R 


-0.038 ± 0.005 


-0.015 ±0.007 a = 0.29 


(B-R)/B 


-0.032 ± 0.005 


-0.007 ±0.007 a = 0.24 


(B-V)/V 


-0.032 ± 0.004 


-0.015 ±0.005 a = 0.35 


(B-V)/B 


-0.027 ± 0.004 


-0.007 ±0.005 a = 0.30 


(V-R)/R 


-0.009 ± 0.002 


-0.004 ±0.002 a = 0.21 


(V-R)/V 


-0.010 ±0.002 


-0.004 ± 0.002 a = 0.24 



